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Spiral waves are investigated in chemical systems whose underlying spatially- 
homogeneous dynamics is governed by a deterministic chaotic attractor. We show 
\£) . how the local periodic behavior in the vicinity of a spiral defect is transformed to 

Q\ ' chaotic dynamics far from the defect. The transformation occurs by a type of period 

0^ . doubling as the distance from the defect increases. The change in character of the 

dynamics is described in terms of the phase space flow on closed curves surrounding 
' the defect. 

Spiral waves are commonly observed in oscillatory and excitable media. |j| They are often re- 
sponsible for the patterns one sees in chemical systems and can give rise to spatiotemporal states 
such as defect-mediated turbulence whose erratic dynamics is characterized by the creation and de- 
struction of pairs of defects (spiral wave cores) with opposite topological charge. M The topological 



^-j- ■ charge n t is defined by 

o 
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V0(r) • dl = th , (1) 



o . 

^sO , where <p(r) is the local phase and the integral is taken along a closed curve surrounding the defect. 

In this article we examine the nature of such spiral wave states in chemical media where the 
underlying dynamics is itself chaotic. More specifically, we consider systems where the dynamics 
^■y of the spatially homogeneous system, described by ordinary differential equations, has a determin- 

istic chaotic attractor which arises through a period-doubling cascade. Consequently, the simplest 
models for the dynamics considered here require at least three phase space variables in contrast to 
the two- variable descriptions of excitable or oscillatory media. jlj We examine how this local de- 
terministic chaos can support spiral waves in the spatially-distributed medium. We also show that 
O ■ systems of this kind can exhibit defect-mediated turbulence and demonstrate that the underlying 

local temporal dynamics is quite different from that in simple oscillatory media. 
Consider the Willamowski-Rossler p| reaction-diffusion equations, 
rS , / 

= KiC x - K-icl - K2C x C y + K-l(? y ~ KiC x C z + K_ 4 + DS/ 2 C x , 
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K2C x C y - K_ 2 Cy - K3C y + K_3 + DV 2 C y , 

-K 4 C X C Z + K_4 + K5C Z ~ K-5C 2 Z + DV 2 C Z , (2) 



where c T (r,t) is the local concentration of species r = x,y, z (we have suppressed the arguments 
of c T on the right hand side of (^|)) and n±i are rate coefficients that contain the concentrations of 
species that are fixed to maintain the system out of equilibrium. The diffusion coefficients of all 
three species are equal to D. 

Suppose the system is spatially homogeneous and the dynamics is described by the ordinary 
differential equations based on the reactive terms in (|J). The resulting mass action rate law 
supports a chaotic attractor that arises by a period-doubling cascade. QJ^] The chaotic attractor 
in the (c x c y c z ) phase space is oriented so that its projection onto the (c x , c y ) plane clearly exhibits 
the (folded) phase space flow around the unstable fixed point (focus) c* = (c* , c*, c*) which spawned 
the attractor. Consequently, to define the phase angle <ft we change variables from c = (c x ,c y ,c z ) to 
a cylindrical coordinate system (p, (j>, z) with origin at c* and z directed along c z . || As the system 
undergoes a sequence of period-doubling bifurcations <p increases by 27r2 ra with each period of the 
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oscillation, where 2 n , n = 1, 2, 3 . . . oo, is the periodicity of the attractor. For period- 1 oscillations 
the other two variables can be uniquely parametrized by the phase <fi, p = p(4>), z = z(cj)) but 
this is no longer true after the first period doubling bifurcation. However, this variable suffices for 
the determination of the location and charge of a topological defect in the spatially-distributed 
medium. 

Next we consider the spatially-distributed chaotic system. Figure [l] is a plot of the local phase 
angle 0(r) in a two-dimensional medium obtained by numerically integrating (Q). Q| One sees a 
complex pattern of spiral defects whose number varies with time. All defects have topological 
charge rit = ±1. The system evolved from an initial state with a single defect in the center of the 
system ||. In the early stages of the evolution the overall number of defects grows rapidly and then, 
depending on the ratio of the diffusion length to the system size, it either saturates monotonously 
or rises to a maximum and decreases to some stationary average value about which it fluctuates. 
Note the different rates of evolution in certain parts of the medium. While the dynamics has an 
almost periodic character in the regions subject to the organizing influence of large, well-established 
spirals (upper left and lower right corners of the panels) the evolution is much faster in domains 
where vortex-antivortex birth and annihilation take place (the vortex-antivortex pair seen in the 
lower left corners of the first two panels disappears as time increases). From visual inspection of 
the individual snapshots of the local phase in this figure it is difficult to detect differences between 
this type of defect-mediated turbulence and that in oscillatory media. Nevertheless, fundamentally 
different kinds of local dynamics consisting of perturbed period-doubled cycles and chaotic motion 
underlie and influence the dynamics of the spiral structures seen in this figure. 

The very fact that stable spiral waves exist in a medium with underlying chaotic dynamics 
demonstrates that the reaction-diffusion kinetics in the vicinity of the spiral centers is by no means 
chaotic - the spiral dynamics locally suppresses the chaos. J^] However, the behavior in those parts 
of the medium which are sufficiently far from any topological defect shows neither temporal nor 
spatial order. 

Thus, we must consider how the local, period-1, spatio-temporal dynamics near a defect is 
transformed into complex chaotic dynamics far from a defect. For this purpose we now examine 
the local dynamics in a reference frame that is centered on one of the defects. We denote the 
time-dependent position of a defect in the xy plane as Td{t) and work in a polar coordinate system 
(r, 8) centered on Yd(t). The local concentration fields may now be expressed in these coordinates: 
c(r,8,t). To investigate this local dynamics in detail we choose the diffusion coefficient to be 
sufficiently large and use no-flux boundary conditions so that the medium supports a single spiral 
wave that persists for long periods of time. Panels (a)-(c) of Fig. || show trajectories in the 
concentration phase space, c(r, 8, t), for several values of r along a ray emanating from the defect 
oriented at an angle 8 — 3/4ir. One observes a period-doubling progression from a perturbed 
period-1 limit cycle near the defect, to perturbed period-2 and period-4 attractors as the distance 
increases. This basic pattern of period doublings is observed for all angles but the circular 
symmetry of the system is not maintained. The origin of this effect will be discussed below when 
the dynamics on closed paths surrounding the defect is examined. 

Variants of this phenomenon were seen in the defect-mediated turbulent regime when the dy- 
namics was observed in a coordinate frame centered on a moving defect. Far enough from the 
defect, locally one finds chaotic dynamics (cf. Fig. || (d)) which is suppressed in the vicinity of 
the defect. Viewed from the defect, as r increases, chaos appears by a truncated period-doubling 
cascade, similar to that seen in stochastically perturbed versions of flows and maps with no spatial 
degrees of freedom, jll] We arc now faced with the question of how the locally periodic dynamics 
can transform to locally chaotic dynamics as r varies and maintain the spiral wave structure in the 
vicinity of the spiral core. 

The nature of this transformation can be deduced from the temporal behavior of c(r, 8, t) for 
some fixed value of r, say r'. Suppose the defect has n t — +1 so that the phase (j> must change by 
2ir as 8 varies through 2ir and r' is such that the local dynamics at (r', 8) is a period-2™ attractor. 
Thus, even though the phase space trajectory of an individual point on the circle with radius r' 
may be a period-2 n attractor, the set of points S — {c(r,8,t) : < 8 < 2tt, r = r' , t = ti] 
must form a closed curve S = S{c) which loops once around c* in the c phase space. The curve 
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S cannot span the period-2™ attractor which loops c* 2™ times in the course of a cycle. However, 
S deforms as time evolves in such a manner that the trajectory of a point on S may be that of a 
period-2" orbit, but the continuity of S is maintained throughout the evolution. 

To make the nature of this deformation of S clear, consider Fig. || which shows S = {c(r, 8, t) : 
< 8 < 2-7T, r = 0.219, i] for t = t 1 ,t 2 ,t 3 ,t4 where r is measured in the fractions of the system 
size. For r — 0.219 and almost all 8 the local attractor is a period-2 limit cycle (P-2) with a small 
inner loop and a large outer loop. For reference, one period of P-2 is shown in each panel of Fig. |^ 
as a light solid line. The filled diamonds are points on S. One sees that at any fixed time instant, 
t = ti, S is a simple closed curve in the c phase space. This curve deforms as time progresses as 
follows: In||(a) S is a large closed curve that lies on the outer loop of P-2. As time increases (b) 
S deforms so that its upper left portion lies on the inner loop of P-2 and smoothly joins to that 
portion of <S remaining on the outer loop. In panel (c) at time £3 one sees that deformation is 
complete and S lies entirely along the inner loop of P-2. Finally, in (d), S expands so that when 
the expansion is complete it again lies along the outer loop of P-2 as in panel (a). Analogous but 
more complicated versions of such deformations are seen for higher periodic orbits. 

We have observed above that in the course of the dynamics the curve S joins portions of the 
inner and outer loops of the P-2 attractor (panel (c) in Fig. |j). In view of the continuity of the 
medium, this implies that there must be a point on a circle with r = r' where the period-2 bands 
merge and the dynamics becomes effectively period- 1. This accounts for the above-mentioned 
broken circular symmetry and the fact that the period-doubling progression observed along rays 
may exhibit different characteristics due to band merging at certain angles 8. 

As one moves from the vicinity of a defect to the bulk medium the nature of the spatiotemporal 
dynamics may be summarized as follows: Close to the defect the local dynamics is approximately 
periodic and as r increases a critical value of r is reached where the local dynamics undergoes a 
"bifurcation" to period-2. This process continues until a chaotic attractor is obtained. As in noisy 
maps and flows this spatial period-doubling sequence truncates at some finite value beyond which 
a noisy chaotic attractor is found. The number of noisy period doublings that may be observed is 
a function of the system parameters, for example the diffusion coefficient, and the characteristic 
distance between defects or between a defect and the boundaries. 

The phenomena described here should exist in any system exhibiting a period-doubling se- 
quence to chaos and be experimentally observable, for example, in chemical reactions carried 
out in continuously-fed unstirred reactors. ||l2fl If conditions are adjusted so that the spatially- 
homogeneous system supports a chaotic attractor, well known in the Belousov-Zhabotinsky re- 
action (l3| , then suitable initial conditions, similar to those that are commonly used to initiate 
spiral waves in excitable media, should produce the spiral wave states described here. Typically, in 
experiments on well-stirred systems, the period doubling sequence is often difficult to resolve since 
it occurs in a very narrow parameter range. This does not imply that the period doubling that 
occurs as one moves away from the defect is confined to a narrow spatial domain. To observe the 
phenomena described in this paper it is only necessary to place the system in the chaotic regime, 
or near it. The spiral structure will locally organize the dynamics and the passage to the chaotic 
dynamics far from the defect should be observable. 

The work suggests the possibility of a variety of phenomena whose existence depends on at least 
three phase space variables. Questions concering the nature of defect dynamics and interactions 
in chaotic media remain to be explored. 
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FIG. 1. Defect-mediated turbulence in a chaotic medium. The local phase <f>(x,y,t) is shown as gray 
shades in the spatial xy plane. Defects can be located as the termini of the white, equiphase contour lines. 
The time interval between frames corresponds to one period of the spiral rotation. Time increases from left 
to right and top to bottom. Rate parameters are: K\ — 31.2, k_i = 0.2, K2 = 1.572, k_2 = 0.1, K3 = 10.8, 
k_3 = 0.12, k 4 = 1.02, k_4 = 0.01, k 5 = 16.5, k_5 = 0.5. The integration time step is At = 5 x 10 -4 and 
the scaled diffusion coefficient is DAt/(Ax) 2 = 10 -3 . Periodic boundary conditions are used. 
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FIG. 2. Phase space trajectories c(r 7 8,t) for K2 = 1.567 and the other rate parameters as in Fig. |l| at 
fixed 6 = 3/4-7T for three values of r: (a) period-1 limit cycle at r = 0.194; (b) period-2 orbit at r = 0.206 
and (c) period-4 orbit at r = 0.344. The scaled diffusion coefficient is DAtJ(Ax) 2 = 0.01 and no-flux 
boundary conditions are used, (d) Local chaotic dynamics at a point in Fig. hi far removed from a defect. 
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FIG. 3. Plots of S in the c phase space for r = 0.219 for four different times tc (a) t\ — 0.0; (b) 
t% — 0.209; (c) t$ = 0.326 and (d) = 0.721; time is expressed in the fractions of average time period 
of P-2. S is shown as filled diamonds while the solid curve represents one period of the P-2 orbit at 
(r = 0.219, 8 = 3/471"). The simulation conditions are the same as in Fig. bl 
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